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Abstract 

We associate a signed digraph with a list of matrices whose dimensions permit 
them to be multiplied, and whose product is square. Cycles in this graph have a 
parity, that is, they are either even (termed e-cycles) or odd (termed o-cycles). 
The absence of e-cycles in the graph is shown to imply that the matrix product 
is a p)-matrix, i.e., all of its principal minors are nonnegative. Conversely, the 
presence of an e-cycle is shown to imply that there exists a list of matrices 
associated with the graph whose product fails to be a Po-matrix. The results 
generalise a number of previous results relating P- and p)-matrices to graphs. 
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1. Introduction and statement of the main results 

P-matrices are square matrices, all of whose principal minors are positive. 
Po-matrices [l[ are square matrices all of whose principal minors are nonneg- 
ative, i.e., matrices in the closure of the P-matrices. We will be interested in 
real Po-matrices. A well-known result of Gale and Nikaido [2|, whose extensions 
and corollaries are discussed in Q, states that a differentiable function on a 
rectangular domain in R™ with P-matrix Jacobian is iniective. This result has 
a number of practical applications - see for example gg. 

The results of Gale and Nikaido have various graph-theoretic implications, 
explored in [f| 0, H| for example. The key results in these papers involve de- 
termining sufficient graph-theoretic conditions for a set of matrices to be Po- 
matrices, followed by additional "nondegeneracy" conditions which guarantee 
that they are in fact P-matrices. Here, a general result is developed, of which 
some of these previous results become corollaries. 

Let k be an arbitrary positive integer. Throughout this paper, a subscript 
or superscript j assumed to belong to {0, . . . , k — 1} should be read as j mod k. 
Let no, ... , rik-i be positive integers. For each j <E {0, . . . , k — 1}, let be an 
i matrix, and define the uq x no matrix A = A^A^ ■ ■ ■ A^ -1 ). We 
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will associate with [A^°\A^, . . ., A^ -1 )] a signed digraph G^(o)...^(fc-i) , which 
will belong to a category of graphs termed signed (k, {1})-BC digraphs, to be 
defined below. The structure of these graphs will imply that all its cycles have 
length which is a multiple of k. 

Given a cycle C with kr\ edges, of which ri have negative sign, we define C 
to be an e-cycle if (— l) r i+ r 2 _ \ anc j an o-cycle otherwise. A signed (k, {1})- 
BC digraph containing no e-cycles will be termed "e-cycle-free" . The first main 
theorem in this paper is: 

Theorem 1. If G = G Am ... A(k -i) is e-cycle-free then A = A^A^ ■ ■ ■ A^'^ 
is a Po-matrix. 

A matrix M determines the qualitative class Q(M) [§] consisting of all ma- 
trices with the same sign pattern as M. Explicitly, Q(M) consists of all ma- 
trices X with the same dimensions as M, and satisfying My > => Xy > 0, 
My < => Xij < and My = => = 0. Given two matrices M and N of 
dimensions such that they can be multiplied, we write 

Q{M)Q(N) = {M'N \M' G Q{M),N g Q{N)}. 

This definition extends naturally to any ordered set of multiplicable matrices. 
Note that Q(M)Q(N) is in general different from Q(MN), and is often not a 
subset of any qualitative class. The second main theorem in this paper is: 

Theorem 2. All matrices in Q(A (QS >)Q(A (lS >) ■ ■ ■ Q{A {k -^) are P -matrices if 
and only if G A {o)j l w---A( k - 1 ) * s e-cycle-free. 



2. Signed (k, {1})-BC digraphs 

Consider a digraph G with vertex set V(G) and edge set E(G). Let S be 
any set of integers. G will be termed "(fe,5)-block circulant", abbreviated to 
(fc,5)-BC, if 

1. V(G) is partitioned into k sets {Vo, Vi, . . . , Vk-i}. 

2. For k > 2, if (j — i) mod k £ S, then there is no (directed) edge from a 
vertex in Vi to a vertex in Vj. 

Every digraph is trivially a (1,5)-BC digraph for arbitrary 5*. When k > 2, 
(fc, 5)-BC digraphs are a generalisation of circulant digraphs ( 13, U| for exam- 



ple). Note, however, that vertices in a (k, S)-BC digraph may have arbitrary 
outdegree and indegree. Here, only the special case S = {1} concerns us. In a 
(fc, {1})-BC digraph, a (directed) path from a vertex in Vj to a vertex in Vj must 
include a vertex from each Vj/ , j ^ j. It follows that all cycles in a (k, {1})-BC 
digraph have length which is a multiple of k. 

Remark. Although, for k > 2, (k, {1})-BC digraphs are fc-colourable, k is 
not in general the chromatic number of G: for example, a (2r, {1})-BC digraph 
with nonempty edge-set is in fact bipartite. 
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Let Vj contain rij vertices. Assume some ordering on these vertices and let 
(i 6 {1, . . . , rij}) refer to the ith vertex in Vj. As usual, an edge (v, v) refers 
to the edge directed from v to v. 

A digraph G is signed if there is a function sign : E(G) — > { — 1,+1}. A 
signed (k, {1})-BC digraph G = G A {t» ... A (k-i) is associated with a list of matri- 
ces [A^> , . . . , A^ fe_1 ^] as above in a very simple way: (i) for each j = 0, . . . , k — 1, 
\Vj\ = rij, (ii) there exists an edge (VJ, Vf +1 ) in G if and only if (A^) rs =/= 0, 
and (hi) the edge (VJ, Vf+i) takes the sign of (A^ J ') rs . Note that entries in the 
matrices are in one-to-one correspondence with edges in G, and that the 
sign-pattern of A^ is in fact a block in the (signed) adjacency matrix of G. 

Example. As an example consider the matrices: 

where a, b, c, d, e, /, g, h, w, x, y and z are arbitrary positive real numbers. Asso- 
ciated with the product A^A^A^ is the signed (3,{1})-BC digraph shown 
in Figure [T] 




Figure 1: The signed (3, {1})-BC digraph corresponding to the product of the three matrices 
J 4(0)^( 1 )^4( 2 ) ; n Eq.[T] The graph has been laid out to emphasise its circulant structure. The 
vertices in Vb,Vi and V2 are represented as filled circles, open circles and boxes respectively. 
Negative edges are represented as dashed lines while positive edges are bold lines. Labels on 
each edge represent the absolute values of the corresponding entries in the matrices, and are 
not strictly part of the graph: they have been added to indicate the correspondence between 
edges and matrix entries. 
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Defining O mxn to be the m x n zero matrix, the signed adjacency matrix of 
the graph G in Figure Q] has block structure: 







2x2 



V 



03x2 

1 1 

-1 1 



-1 



^3x3 
^2x3 







2x2 




It can be seen immediately that each block is simply the sign-pattern of A^°\ 

Although the graph in Figured] has a number of cycles (both of length 3 and 
of length 6), all of these can be computed to be o-cycles, and so, by Theorem [TJ 
the product A^A^A^ 

is a Po-matrix. This is true whatever the magnitudes of 
the entries in the matrices. Clearly, given isomorphic signed digraphs G\ = G2, 
Gi is e-cycle-free if and only if G2 is e-cycle-free. Since 

G A tl) A (2) A (o) = G A (2) A (0) A (l) = G A (0) A (l) A (2) , 

^(1)^(2)^4(0) and ^4(2)^4(0)^(1) are also p _ matrices . 



3. Preliminaries needed for the proofs 

Permutations. Consider an ordered set a — [eti, a>2, ■ ■ ■ , ot r ], and a bijec- 
tion (3 : a — > a. Defining fa = j3(ai), the ordered set [fa, 02, ■ ■ ■ , fa] will also be 
referred to as fa It will always be clear from context whether an object referred 
to is a bijection or an ordered set. 

Any permutation j3 has a parity P((3), i.e., P{fa) — +1 if (3 is an even 
permutation and P(/3) = —1 otherwise. Given two permutations /? and /3 , 
P(fa3') = P(/3)P(/3'), implying that P{fi) = P(/3" 1 ). Note the following ele- 
mentary result about the parity of permutations. 

Lemma 3. Consider a permutation (3 of a finite set of size r. Write P as a 
product of disjoint cycles, C\, . . . , C s (1 < s < r), including trivial cycles. Then 

p(fa = (-iy-> 

i.e., [3 is even (resp. odd) if the total number of elements in (3 minus the total 
number of cycles in its decomposition is even (resp. odd). 

Proof. See [Hj], for example. 

From here on a« will always refer to a nonempty subset of {1, . . . ,nj}, 
and will be assumed to have the natural ordering, am will refer to the mth 
element in so that < c4 < a 3 < ' ' ' • Given some , define 
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= \yk | g a^')} C Vj. will refer to a permutation of a^. Given the 
one-to-one correspondence between the elements in , and vertices in V J ai3) , 
can equally be regarded as a permutation on IA Q<J> . 
Now consider some sequence (a^, . . . , cr fe_1 )), such that \a^\ = \a^'\ 
for each and a corresponding sequence of permutations ((3^°\ ■ ■ . ,/3^ fe_1 ^). 
Define i? = [1, 2, . . . , |a^|], and define the bijections ij : V°' U} — > R by 

m ) = m. In other words, ij associates with each vertex in its order. 

fi^i) then induces the bijection : R — > R defined by fi^i) = lj o (3^ o ij 1 . 
Further, define the bijections <j)j : V" 3 — > V^ +1> by 4>j = lJ^ x o o ij, i.e. 

^jv^j m ) — ^j+i ■ Equivalcntly, cpj = fj^ +L > o o t,-. Figure [2] illustrates 
all of these relationships. Note that in the special case k = 1, = 




Figure 2: The commutative diagram which encapsulates the relationships /8") = OjofiU) oi 
fa = /3O+ 1 ) o iT^ o tj, and ^ = iT^ o /3O+ 1 ) o Lj . 



4. Proof of Theorem Q] and its immediate consequences 

The following notation is used. Given an r x s matrix M , and two (nonempty) 
ordered sets 7 C {1, . . . , r} and S C {1, . . . , s}, then M(7|<5) is the submatrix 
of M with rows indexed by 7 and columns indexed by d. If I7I = \S\, then 
M[y\S\ = det(M(7|<5)). We write M[y] as shorthand for M[y\^}. If 7 and S are 
of equal size, then M 7) $ will refer to ^y.A- 

PROOF of Theorem [1] The cases k = 1 and fc > 2 are conceptually similar: 
however in order to avoid notational difficulties, they are presented separately. 

Case 1: k = 1. Let V = V(G). Choose and fix some nonempty a C 
{1, . . . , n}, and consider the minor A[a\. Given a permutation /3 of a, define T, 
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a term in A[a], by: 

T = P(/3)A a ,p . 

Suppose that T ^ 0. Since nonzero entries in A are in one-to-one correspon- 
dence with edges in G, corresponding to T is an edge-set E C E(G) defined 
as follows: the edge (V r ,V s ) is in E if and only if (r, s) = (a m ,/3 m ) for some 
integer m. Equivalently, the edge (V r , V s ) is in E if and only if V s — f3(V r ). 

The set of endpoints of edges in E is precisely V a . Consider the ver- 
tex V° m . Then the incoming edge (/3~ 1 (V rQm ), V° m ) and the outgoing edge 
{V a ™ , /3(V a ™ )) are the only edges from E incident on V a ™ . If (3{V am ) = V a ™ , 
then these edges coincide and in fact there is a loop at V am . Otherwise, the 
edges are distinct. In either case, E, regarded as a subgraph of G, consists of 
vertex- and edge-disjoint cycles. 

Suppose p m (v) = v for some vertex v, but j3 q {v) ^ v for q < m. Then 
the vertex v lies on a cycle in E of length m. Thus, decomposing /3 as a 
product of disjoint cycles (including trivial cycles), these cycles are in one-to- 
one correspondence with cycles - in the graph-theoretic sense - in E. Trivial 
cycles correspond to loops. Assume that there are N such cycles. By LemmaO 
P(0) = (-l)M-tf 

Returning to the term T in the expansion of A[a], 

sign(T) = P(/3)sign04 Q)/3 ) (2) 

Consider a cycle C in E including n edges of which ri are negative, and define 
sign(C) = (— l) T2 (i.e. sign(c) is the product of signs of edges in C). If C is 
an e-cycle, then (— l) ri+r2 = 1, and so sign(C) = (— l)' 1 . Similarly if C is 
an o-cycle, then sign(C) = (— l) ri+1 . Decompose E into disjoint cycles, which 
comprise N e e-cycles, and N a o-cycles (so that N = N + N e ). Let 8 be the 
total number of edges in o-cycles and 8 e the total number of edges in e-cycles, 
so that 8 + 8 e = \a\ (since there are \a\ edges in E). Taking the product of 
signs of edges in E over e-cycles and o-cycles separately gives 

sign(A^) - (-l) e °(-l)°'+ N ° - . (3) 

Since P(/3) = (-1)1°!-^, and sign (A Qj(9 ) = (-1)H +A, °, Eq. Ogives: 

sign(T) = (_i)l«l-A r (_i)l«l+Wo = (_!)*.. (4) 

If G is e-cycle-free, then N e = in this expression, in which case sign(T) = 1. 
Since T is an arbitrary nonzero term in A[a], A[a] > 0. Since a is an arbitrary 
nonempty subset of {1, . . . , n}, A is a Po -ma trix. 

Case 2: k > 2. Choose and fix some nonempty C {1, . . . , no}, and 
consider the minor ^4[a^ ^]. Enumerate all lists (a*- 1 -*, . . . , a^" 1 ^) where for each 
i = 1, . . . ,k — 1, aW i s a nonempty subset of {1, . . . , rij} such that \a^ 1 ' \ — \a^°'\. 
There may of course be no such subset for some i, and hence no such lists. 

By the Cauchy-Binet formula [3] applied recursively, 
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The sum is over all possible lists (a^- 1 a^ k 1 ' > ) with \eft'\ — \a^°'\. Now 
choose and fix some particular choice aS 1 ', . . . , a'-' " 1 -', and choose permutations 
fi \. . .,^ k ~ 1 \ For each j define T h a term in A^[a^\a^+% by: 

Suppose that for each j, Tj is nonzero so that T = J| ■ Tj ^ 0. Note that T 
is then a nonzero term in the expansion of A[a^ ^]. Since nonzero entries in the 
matrices are in one-to-one correspondence with edges in G, corresponding 
to T is an edge-set E C E(G) defined as follows: the edge (VT, Vf +1 ) is in E 

if and only if (r, s) = (a$,j8& +1) ) for some integer m. Equivalently, the edge 
(VJ, V? +1 ) is in E if and only if (j>j(VT) = Vf +1 . 

The set of endpoints of edges in E is precisely [Jj , and in fact each 
such vertex has exactly two edges from E incident on it, one incoming and 
one outgoing. For example, consider the vertex VJ", where r S ofc>>. Then the 
incoming edge ((f)J_ 1 (VJ') , Vf) and the outgoing edge (Vf , <j)j(V r )) are distinct 
edges in E, and are, by the definition of E, the only two edges in E incident on 
Vj". As a consequence, E, regarded as a subgraph of G, consists of vertex- and 
edge-disjoint cycles. 

Next, consider the bijection 4> : V^"' — > Vq {0) defined by = 4>k-i ° 0fc-2 ° 
• • • o (p . Suppose 4> m (v) — v for some vertex v, but <f> q {v) ^ v for q < m. 
Then the vertex v lies on a cycle in E of length km. Decomposing </> as a 
product of disjoint cycles (including trivial cycles), these cycles are in one-to- 
one correspondence with cycles - in the graph-theoretic sense - in E. Assume 
that there are N such cycles. By Lemma|3l P{4>) — (— 1)I Q<0> \~ N . Applying the 
relations 4>j — L J+i ° /3^ +1 - ) o ij gives 

so that = n s - n -ro 1 

Returning to the term T in the expansion of ^a 1 - - 1 ], 

sign(T) = nsign(T ? -)= ( JJ P^ 1 )) I ( JJ sign ] . (5) 

J=0 \j=0 / \i=0 / 

The first term in this expression has already been determined: from above, 

n^°' +i) ) = f[pw {j) ) = p{<t>) = (-l)"^"-^ . (6) 

3=0 3=0 
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Consider a cycle C in E including kr\ edges of which r-i are negative. As in the 
case k = 1, if C is an e-cycle, then sign(C) = (— l) n , while if C is an o-cycle, 
then sign(C) = (— Decompose E into disjoint cycles, which comprise N e 
e-cycles, and N a o-cycles. Let kOo be the total number of edges in o-cycles and 
k9 e the total number of edges in e-cycles, so that O + 9 e = \a^\ (since there 
are fc|cr°)| edges in E). Taking the product of signs of edges in E over e-cycles 
and o-cycles separately gives 

n W = = (-i)'« (o, ' + ^ . ( 7 ) 

3=0 

Substituting Eqs. [5] and [7] into Eq. [3] gives: 

sign(T) = (_l)l« (0) l-"(_l)l« (0) l+"o = (_!)*.. (8) 

Note that this is just Eq. 2] again. As in the case k = 1, if G is e-cycle-free, then 
sign(T) = 1, and since T is an arbitrary nonzero term in ^[a' - 1 ], ^[o^ * 1 ] > 0. 
Since a*- - 1 is an arbitrary nonempty subset of {1, ...,no}, A is a Po _ma trix. 

□ 

We have the following corollary to Theorem [T] 

Corollary 4. Consider a square matrix A = A^A^ ■ ■ ■ A {k -^ such that G = 
G A m A (i)... A{k -i) is e-cycle-free. Then matrices in Q(A^)Q(A^) ■ ■ ■ Q{A^ k -^) 
are all P^-matrices. 

Proof. By definition, any matrix B E Q(A (0) )Q(A (1) ) • • • Q(v4 (fe_1) ) can be 
written B '= B^B^ ■ ■ -P^ 1 ), where G Q(A«). But 

Gb(o)b(i)..._b(*-i) — GaWAW—AI*- 1 )) 

so, by Theorem [TJ B is a Po-matrix. □ 
Remark. For each r = 1, . . . , k — 1 the matrix product 

A (r) A (r+l) , , . ^(fe-l)^(O) . . . ^(r-l) 

gives rise to a graph isomorphic to G = G A (o) A (i) ... A (k-i) . Thus, in fact, if G is 
e-cycle-free, then all matrices in 

Q{A^)Q(A^ r+i y) ■ ■ ■ Q(A^)Q(A^) ■ ■ • Q(A (r_1) ) 

are Po-matrices. 

5. Converse results and proof of Theorem [2] 

A variety of converse results are possible, that is, results which guarantee 
that if a signed (k, {1})-BC graph contains e-cycles, then there exist matrices 
in some set which fail to be Po-matrices. The most useful formulations depend 
on the application. Lemma O is the basic result from which such results follow: 
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Lemma 5. Consider a graph G — G y i(o)^(i)... J 4(fe-i) such that all edges in G 
lie on a single e-cycle C of length kr. Then A = A^°'A^ • ■ ■ A^ k ~^ is not a 
Pq- matrix. 

Proof. Define the sets (a^°\ . . . , a^ k ~ x ">) by the stipulation that s G aW if and 
only if Vf lies on C. Clearly \a^\ — r for each i. Each vertex in V™ 3 lies on 
C and hence has exactly two edges, one incoming, and one outgoing, incident 
on it. So it is possible to define bijections <pj : V" — s- V-* +1 as follows: given 
vertices v E V" 3 and v <E Vj+^ % , 4>j( v ) = v if there is a directed edge (v, v) 
in E. Each </>j induces a permutation : V?j_i +I — > V?j_i +1 defined by 

^O'+i) = ^ o tr 1 o (see Figure [2]). 

Consider the minor A[a { % There is a nonzero term in this minor 

j=0 3=0 

Moreover T is the unique nonzero term in j4[o/°)] : any other nonzero term would 
imply the existence of an index j and a permutation S : V°^i ' — > ' , 

S ± such that A^ u) s ^ 0. Letting s be an index such that S s ^ P { s +1 \ 

A^y) must then be nonzero, implying the existence of an edge {V^ s , Vj+i) 

in G which does not lie in C. But by assumption C contains all edges in G. 

By Eq. [5J sign(T) = (— 1) e , where N e is the number of e-cycles in the 
subgraph associated with T. Since this subgraph is precisely C, N e — 1 and 
sign(T) = -1. Thus A[a (0) ] < and A fails to be a P -matrix. □ 

Corollary [S] illustrates an application of Lemma [SJ 

Corollary 6. Suppose a graph G^(o)yt(i)...^(<j-i) contains an e-cycle C. Then 
there are matrices in X = Q{A^)Q(A^) ■ ■ ■ Q(A ( - k ~ 1 ^) which are not Pa- 
matrices. 

Proof. Each edge in C corresponds to an entry in one of the matrices A^\ 
For each j define A^> to be the matrix A^ with all entries not corresponding to 
edges in C set to be zero. Then the matrix factorisation A = A( MW • ■ • 
gives rise to a graph which consists solely of the e-cycle C, and hence, by 
Lemma A fails to be a Po- ma trix- But A G cl(A') (that is the closure of 
X), and since the set of Po- ma trices is closed, there are matrices in X which fail 
to be P . □ 

Proof of Theorem [2j This is immediate: Theorem[2]is simply a combination 
of Corollaries Hand H □ 
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6. Notes and conclusions 



Although the special case k = 1 was treated for completeness, the result can 
easily be inferred from previous work. For k — 1, Theorem [5] states that given a 
square matrix A, all matrices in Q(A) are Po-matrices if and only if the (unique) 
signed (1,{1})-BC digraph Ga associated with A is e-cycle-free. However Ga 
is closely related to an object often called the interaction graph or I-graph in 
the literature. In fact the I-graph associated with A is just Gat- Results in 
[H, [|| showed that G_a t lacks positive cycles if and only if all matrices in Q(A) 
are Po-matrices. Trivially, G_^r lacks positive cycles if and only if G-a lacks 
positive cycles. The definitions imply that e-cycles (resp. o-cycles) in Ga are in 
one-to-one correspondence with positive cycles (resp. negative cycles) in G-a- 
So G-a lacks positive cycles if and only if Ga is e-cycle-free. Together these 
observations imply that matrices Q(A) are all Po-matrices if and only if Ga is 
e-cycle-free. 

The case k = 2 has also effectively been treated previously in 0, Q , where the 
associated graphs were termed "DSR graphs". The main differences between 
the definition of a DSR graph in [8|, and a signed (2, {1})-BC digraph here, are 
(i) directions on all edges are reversed, (ii) here, edge-labels have been ignored, 
while some computations on DSR graphs in [3, |8[ involved edge-labels, and (iii) 
in the construction of the DSR graph a pair of identically signed edges (v,v) 
and (v,v) are replaced with a single undirected edge, with a view to removing 
o-cycles of length 2 from the graph, and thus simplifying computation. This 
process neither creates nor destroys e-cycles, and so does not change the key 
fact that an absence of e-cycles implies that associated matrices are Po-matrices. 

The treatment in @,@] also suggests that extensions obtaining sharper results 
by introducing edge-labelling and more complex computations on the graphs are 
possible. The most useful forms that such extensions might take depend on the 
applications in question. These directions will be treated in future work. 
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